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tion on the Galactic and extragalactic interstellar medium along 
the lines of sight to the two neighbouring galaxies. 

A first analysis of the ISM in the direction of SN 1993 J us- 



ing lUE spectra was made by de Boer et al. (1993). Optical 



1994). High velocity clouds (HVCs) and intermediate velocity 



Abstract. We present an extended analysis of the high dis- 
persion lUE spectra of SN 1993 J in M81, based on the im- 
proved data in the lUE final archive. Eight interstellar absorp- 
tion components can be identified in the velocity range from 

-200 kms-i to +200 kms'i. Column densities are deter- spectra of the supernova were'analysed by Vladil^al. ( [1993 , 

mined for Fe II, Mg I, Mg II, O I, S II, Si II, and Zn II. From the 

ion ratios, we estimate electron densities and temperatures of 

the gas in these components. In addition, we rederive column 

density profiles for the higher ionisation stages Al III, Si IV, and 

C IV for the entire Une of sight. 

Of special interest are the absorption components at pos- 
itive velocities, which are associated with the intergalactic 
medium (IGM) in the foreground of M 8 1 . The medium has 
many gas components, as known from radio measurements, 
probably released in dynamical interactions by galaxies in the 
M 8 1 group. The lUE spectra contain the only UV absorption 
information on such gas. We find neutral and ionised gas ab- 
sorption at the same velocity, having a metal content similar to 
a galactic disk. 

We combine the spectra towards SN 1993J in M81 with 
those towards SN 1987 A in the LMC, the directions being es- 
sentially antipodes on the celestical sphere. Comparing the to- 
tal profile with QSO spectra, we conclude that most absorption 
in QSO spectra occurs in diffuse clouds on those sight lines and 
is not similar to absorption by galactic disk gas. However, one 
absorption component towards SN 1993J is clearly very similar 
to some of the absorption components in QSO spectra. 

Key words: supernovae: individual: SN 1993J - Galaxy: halo 
- intergalactic medium - Galaxies: ISM - quasars: absorption 
fines - Ultraviolet: ISM 



1. Introduction 

In recent years two supernovae occurred in nearby galaxies, 
SN 1987 A in the LMC and SN 1993 J in M 81. High resolution 
UV spectra of both supernovae have been taken with the lUE 
satellite in the full wavelength range from 1150 up to 3200A. 
A large number of absorption lines from various metals is lo- 
cated in this part of the spectrum, providing valuable informa- 
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clouds (IVCs) are present. These designations refer to gas with 
velocity not compatible with what is expected from Galactic ro- 
tation. Rather arbitrarily, HVC is used for |urad| > 100 kms^^ 
and IVC for 100 kms^^ > |?;rad| k. 25 kms"^ We extend the 
previous analysis of lUE spectra by using the rereduced spectra 
from the final archive (lUEFAJJ) to determine column densities 
for the metals in the neutral components. We rederive the col- 
umn density profiles for the high ion stages of Al III, Si IV and 
Civ. 

The absorption components at velocities > 50 kms^^ are 
of major interest here. Their velocities are uncommon for the 
Galactic sky in this direction and the absorption most likely 
originates in or near the M 8 1 gala xy gro up (distance of M 8 1 
3.63 ± 0.34 Mpc, Freedman et al. 1994 ). The column densi- 
ties derived are used to estimate the density and temperature in 
these clouds and to find support for the proposed location. Ap- 
pleton et al. ( 198 1 ) give a list of velocities for members of the 
M81 group. M81 itself has a systemic velocity of —34 km s^^, 
with a rotational velocity of —130 kms^^ at the position of 
SN 1993J (Rots [975]). The velocity of M82 is +183 kms^^; 
that of NGC 3077 +19 km s^^. Intergalactic molecular gas was 
found in CO emission southeast of the M 8 1 disk at a velocity 
of —35 km s^^ (Brouillet et al. 1992| ). Also, intergalactic gas at 
velocities of +120 to +140 kms^ is detected in a tidal tail 
region connecti ng M 8 1 and M 8 2, no rth of the M 8 1 center 
(Appleton et al. |l98lt Yun et al. |l993| ). The gas found in the 
SN 1993J absorption spectra fits into this range of velocities. 

Combining the spectra towards SN 1993 J with those to- 
wards SN 1987 A (essentially on the opposite side of the sky) 
allows one to create the absorption spectrum through 3 galax- 
ies. This total spectrum is compared with absorption as seen 
on QSO lines of sight. This provides a possibility to judge the 
origin of the absorptions. 
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Table 1. lUE high dispersion spectra of SN 1993J 



Image 



1993 Date Time Exp. time [min] 



LWP 25238 30 March 06:13 UT 80 

SWP 47386 30 March 07:58 UT 320 

SWP 47394 31 March 05:43 UT 655 

LWP 25250 01 April 03:41 UT 240 



The data used are presented in Sect. |[ The results from the 
analysis of the SN 1993 J data are given in Sect. for the neutral 
gas and in Sect. for the ionised gas, followed in Sect, p^ by a 
discussion of the structure of the SN 1993 J line of sight. 

In Sect. g| we present the combined spectra of the full line of 
sight, followed by the comparison with QSO spectra in Sect. 0. 
Finally, in Sect. |], we summarise our understanding of the na- 
ture of the intergalactic clouds towards SN 1993 J in relation to 
the dynamics of the M 8 1 group of galaxies. 



2. The data 

We analysed the lUE spectra provided by the ESA lUE Fi- 
nal Archive at VILSPA, Spain. Four spectra of SN 1993 J were 
obtained: two each in the short and long wavelength range 
(SWP 47386, SWP 47394, LWP 25238, LWP 25250). The res- 
olution of the spectra is about 20 — 25 km s ^ ^ . The dates of the 
spectra are listed in Table |l]. For a representative selection of 
absorption line spectra see Fig. |l]. Due to the rapid decrease in 
UV flux in the supernova, the net flux is significantly weaker 
for those spectra which were taken later. This leads to a lower 
signal-to-noise ratio, even though these spectra have larger ex- 
posure times. 

The SN 1993 J data from lUE can be combined with pub- 
lished data from other wavelength regions. In particul ar, we 
use sodium column densities derived by Vladilo et al. (1993, 



1994 ) and by Bowen et al. ( 1994 ) from optical spectra. Bowen 



et al. ( 1994 ) also derive column densities for Mg II from HST 
GHRS spectra, at a higher resolution than that o f the lUE. The 
Effelsberg H 1 21 cm spectra (Vladilo et al. 1994 1 are also taken 
into account. 

A much larger number of spectra was tak en for SN 1987 A. 
Referring to the spectral atlas of Blades et al. ( 1988 ), we use the 
short wavelength spectrum SWP 30381, which has unsaturated 
net flux of FN < 40 000 in the echefle orders 78-115. The spec- 
trum was taken on 25 February 1987 with an exposure time of 
30 min. No long wavelength spectrum of SN 1987 A is used. 

The spectra show the typical velocity profile seen along al- 



most all sight lines towards the LMC (Savage & de Boer 1981 



de Boer & Savage 1980): local Galactic gas around kms^ , 
a Galactic IVC and HVC at +60 and +125 kms"\ respec- 
tively, one absorption component of as yet unassigned location 
in front of the LMC at +160 km s"^, and broad LMC absorp- 
tion between +200 and +330 kms^^. 




-300 -200 -100 



Vlsr [km s' ] 



Fig. 1. Interstellar absorption profiles of some transitions from 
the earher spectra SWP 47386 and LWP 25238. Zn ii (2026 A) 
blends with Mg I at positive velocities. The relative error anjn 
in the continuum is in the vicinity of 10%. The "R" denotes a 
reseau mark 

Metal column densities for SN 1987 A were recently de- 
rived by Welty et al. ( 1999 ) using profile fitting. For the h igh 
ionisation states we use the analysis by Savage et al. ( 1989| ). 

The SN 1993 J spectra discussed here are reduced by the 
NEWSIPS algorithms, in contrast to the analysis published by 
de Boer et al. (1993). This provides an improved background 
correction, especially in the higher, closely spaced echelle or- 
ders, an increased signal-to-noise ratio, and a higher spectral 
resolution (Garhart et al. 1997). We smoothed the data using 



a 3 pixel wide boxcar algorithm. In the NEWSIPS spectra the 
zero intensity level of strong absorption lines appears to lie a 
few percent too low. We have generally lifted the zero level 
accordingly. 
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3. Lower ionisation stages in the SN 1993J spectra 

3.1. Metal and H I column densities 

The wavelength scale of the lUE spectra might be affected by 
pointing errors, resulting in a velocity shift for the absorption 
components. We therefore fix the velocity scale by aligning, for 
each transition, the absorption at '^ and ^ —50 kms^^ with 
the corresponding clearly visible emission in the Vladilo et al. 



( 1994) H I 21 cm spectra. All velocities are given with respect 



to the LSR. 

For the identification of the absorption components we use 
the well exposed transitions of Si II at 1304 A and of Fe II at 
2382 A. Due to the improved data reduction, we are able to 
resolve eight absorption components (see Fig. M with the ve- 
locities listed in Table 0. An additional component might be 
present at —170 kms^^ in some transitions, but it is too weak 
for analysis. No absorption components are seen at velocities 
lower than —180 kms^^ or larger than +200 kms^^. Due to 
the low signal-to-noise ratio of the spectra, the centres of the 
absorption components vary slightly and we refer to the com- 
ponents by rounded radial velocities. 

The component D around kms^^ in many cases splits 
into a strong component D^ and a weak D^. These compo- 
nents, however, cannot be clearly distinguished in all transi- 
tions. Components A and B near —130 and —90 kms^^ are 
resolved, but still seem to blend in some transitions. We anal- 
yse them separately now, in contrast to the previous publica- 
tion. The same holds for F and the weaker component G. The 
blending leads to larger uncertainties in the resulting column 
densities, especially for the weaker of these components. Com- 
ponent E was disregarded in the previous publication, but is 
clearly visible in most transitions. 

The equivalent widths of the absorption lines are measured 
using a trapezium fit. Due to the faint continuum no better re- 
sults would be achievable with Gaussian fitting. We account 
for uncertainties in the continuum by determining equivalent 
widths for a lowest and highest possible local continuum. Half 
of this variation is quadratically added to the error resulting 
from noise. Metal column densities are then determined using 
standard curve-of-growth methods. 

We use the transitions of Fell (1608, 2260, 2344, 2374, 
2382, 2586, and 2600 A), Si II (1260, 1304, 1526, and 1808 A), 
Mgl (2852 A), Mgll (2796 and 2803 A), O I (1302 A), S II 
(1250, 1253, and 1259 A), and of Znll (2026 A); / values 



are taken from the compilation by Morton (1991), plus newer 



values from the review by Savage & Sembach (1996). The b 



value for each absorption component is fixed using the equiv- 
alent widths of the Si II and Fe II transitions. The remaining 
equivalent widths are then fitted to the appropriate theoretical 
curve of growth. Column densities and 6-values are given in 
Table g. The errors listed represent minimum and maximum 
column densities possible within the range in b, allowed by the 
errors in the equivalent widths. All subsequent error values fol- 
low from standard error propagation. 

Hi column densities from Hi 21 cm emission spectra 



(Vladilo et al. |1994| ) can only be derived for the components 



C at —50 kms^^ and Da and D^ around kms^^. The broad 
emission from the disk of M 8 1 makes it impossible to give 
reliable Hi column densities for the absorption components 
at higher negative velocities. At positive velocities only upper 
limits follow from the 21 cm spectra. 

However, we can derive H I column densities for all absorp- 
tion components from the column densities of the normally 
only slightly depleted sulfur and zinc, under the assumption 
of no depletion at all for these two metals. Although this as- 
sumption introduces an additional source of error, this method 
provides more reliable information on H I in the single com- 
ponents than the 21 cm radio spectra at large cloud distances, 
since these are based on integration over the full width of the 
radio beam. As shown in Table |[ the H I column densities cal- 
culated this way are in good agreement with the 21 cm mea- 
surements for those components, for which values are given 
by Vladilo et al. (1994). For the other components the column 
densities from the lUE data are equal to or larger than those 
from the radio emission spectra, which is expected for clouds 
at a larger distance that do not fill the radio beam. Based on 
these data, the metal abundances also listed in Table g can be 
calculated. The large depletion of metals in component Df, is 
possibly the result of underestimated metal column densities 
due to blending with the strong component Dq. 

3.2. Physical conditions in the gas 

Important information on the physical conditions in the inter- 
stellar clouds, such as temperatures and electron densities, can 
be inferred from the column density ratios of the atom's differ- 
ent ionisation stages. The various spectra of SN 1993J provide 
column densities for two of these: 

The ratio of iV(Mg ii)/iV(Mg i) can be read directly from 
the lUE spectra for most absorption components. If Mg I col- 
umn densities cannot be determin ed fro m the lUE spectra, we 
take the values of Bowen et al. (1994) from GHRS spectra. 
Since these data have a higher resolution, we have to identify 
their components with those found by us. We assume (veloci- 
ties in kms^^) 

[-145, -115]^ A; [-110, -70] ^B; 

[-70,-30]^C; [-20, +5]^ Da; [+5,+15]^Db; 

[+80, +105] ^E; [+110, +150] ^F; [+150, +175] -» G. 

The resulting values are listed in Table ^ For Mg I in Da we 
assume the larger column density from Bowen et al. (1994) as 
an upper limit. The larger value in E can be due to an additional 
component at 104 kms~^ in the GHRS data. We will use the 
column density from lUE in this component. The other values 
agree within the ranges of the errors. 

Comparing the Mgll profile at 2803 A from lUE and 
GHRS spectra, one finds in the lUE profile what looks like a 
feature between —50 and kms^^ (see Fig. Il|). Considering 
the superior quality of the GHRS spectra, this feature is prob- 
ably not real. The effect of the feature on the resulting column 
densities, however, is negligible. 
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Table 2. Metal column densities log A^ [cm"^] towards SN 1993J 





A 


B 


C 


D, 


Db 


E 


F 


G 


ULSR[knis ^] 


-130 


-90 


-50 


-10 


±10 


±80 


±130 


±160 


Fell 


15.25 ±0.25 


13.90 ± 0.30 


15.15 ±0.20 


14.95 ±0.25 


14.35 ± 0.40 


14.15 ±0.25 


15.30 ±0.20 


14.15 ±0.25 


Mgl 


- 


- 


12.65 ±0.15 


12.45 ±0.20 


- 


11.20 ±0.35 


12.95 ±0.35 


11.80 ±0.15 


Mgll 


14.95 ± 0.30 


> 13.90 


14.15 ±0.15 


14.60 ± 0.30 


< 13.40 


14.70 ± 0.30 


14.60 ± 0.30 


14.20 ±0.30 


Oi 


14.90 ± 0.30 


> 15.85 


- 


- 


13.85 ±0.45 


15.70 ± 0.50 


16.50 ±0.20 


15.25 ±0.45 


Sii 


16.10 ±0.40 


14.60 ± 0.20 


15.55 ±0.25 


15.65 ±0.25 


15.85 ±0.45 


15.30 ± 0.40 


14.25 ± 0.25 


< 14.15 


Sin 


15.05 ±0.15 


14.35 ±0.10 


15.20 ±0.10 


15.05 ±0.15 


14.45 ±0.15 


14.20 ±0.20 


14.10 ±0.20 


14.25 ±0.20 


Znll 


12.90 ± 0.20 


12.35 ±0.25 


12.05 ±0.30 


12.95 ±0.20 


- 


11.85 ±0.25 


12.40 ± 0.20 


- 


blkms-^]" 


6-7 


4-6 


9 


6-7 


1.5-2 


5-7 


6-7 


5-6 



" &- values are derived from curves of growth. Their uncertainties can be up to 2 km s ^, depending on the scatter of the equivalent widths 

Table 3. Hydrogen column densities log A'h [cm^^] and metal abundances log (A^xZ-^Vh) derived towards SN 1993J (discussion 
in Sect. |]) 





A 


B 


C 


Da 


Db 


E 


F 


G 


solar" 


log Nus z„ 


20.60 ± 0.45 


19.58 ±0.32 


19.90 ± 0.39 


20.40 ± 0.32 


< 21.10 


19.78 ±0.47 


19.43 ±0.32 


< 18.95 




logiVH.i™ 




- 


19.81 


20.36 


< 19.4 


< 19.4 


< 19.4 


< 19.4 




Fe 


-5.35 ±0.51 


-5.68 ± 0.44 


-4.75 ±0.44 


-5.45 ±0.41 


> -7.15 


-5.63 ±0.53 


-4.13 ±0.38 


> -5.05 


-4.50 


Mg 


-5.65 ±0.54 


> -5.68 


-5.74 ± 0.42 


-5.80 ± 0.44 


> -7.70 


-5.08 ±0.56 


-4.82 ± 0.44 


> -5.05 


-4.46 


O 


-5.70 ±0.54 


> -3.73 


- 


- 


> -7.70 


-4.08 ± 0.69 


-2, 93 ±0.38 


> -4.15 


-3.16 


Si 


-5.55 ±0.47 


-5.23 ±0.34 


-4.70 ± 0.40 


-5.35 ±0.35 


> -6.80 


-5.58 ±0.51 


-5.33 ±0.38 


> -4.90 


-4.45 


" solar abundances are taken from review by de Boer et al. (|l987a|) 












'' 21 cm H I column densities are taken from Vladilo et al. (199' 
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Table 4. iV(Mgi)/A^(Mgii) ratios 



Table 5. N{Zn ii)/A^(Na i) ratios, based either on the observed 
Zn II and Na I column densities or on Na I alone, using Eq. dl]) 



log iVMg 



log A^Mg I 

(Bowen) 



log A'"Mg li log iVMg i/Nms 



log A^Zn II log iVNa 



logATznll/A^'Na 



A 
B 


12.65 ±0.15 


12.41 ±0.02 
12.00 ±0.08 
12.64 ±0.06 


14.95 ±0.30 

> 13.90 
14.15 ±0.15 


-2.54 ±0.30 

< -1.82 
-1.51 ±0.16 






(Vladilo) 


(Bowen) 






C 


A 


12.90 ±0.20 


12.06 


12.03 ±0.01 


0.86 ±0.20 


obs. 


D, 


12.45 ±0.20 


ismt^ol 


14.60 ± 0.30 


< -0.55 


B 


12.35 ±0.25 


- 


< 11.48 


> 0.62 


obs. 


Di, 


- 


11.41 ±0.03 


< 13.40 


> -2.02 


C 


12.05 ± 0.30 


11.23 


11.34 ±0.02 


0.77 ±0.30 


obs. 


E 


11.20 ±0.35 


11.92 ±0.08 


14.70 ± 0.30 


-3.50 ±0.46 


D„ 


12.95 ±0.20 


12.62 


14.09 ±0.13 


< 0.53 


obs. 


F 


12.95 ±0.35 


12.76«:1° 


14.60 ± 0.30 


-1.84 ±0.33 


Dft 


- 


12.18 


< 10.23 


> -0.84 


Eq.(l) 


G 


11.80 ±0.15 


11.93 ±0.03 


14.20 ±0.30 


-2.27 ±0.30 


E 
F 


11.85 ±0.25 
12.40 ±0.20 


12.89 


< 10.49 
13.01 ±0.01 


> 1.11 
-0.55 ±0.20 


obs. 












obs. 












G 


- 


- 


< 10.32 


> 0.74 


Eq.(§ 



The ratio A^(Naii)/A^(Nai) is not available from observa- 
tions. We can, however, use A^(Znii)/A^(Nai) as a substitute 



(Vladilo & Centurion 1994). The only slight depletion of Znll 



and its ionisation potential of 18.0 eV make it a good tracer of 
the neutral gas, while Na I is sensitive to the physical conditions 
in the gas. 



Column densities for Na I are given by Vladilo et al. (1994 ) 



and by Bowen et al. (1994). We equate the higher resolved op- 



tical velocity components with our lUE components as defined 
above for Mg I. Both optical data sets are quite similar Bowen 



et al. ( |1994| ), however, have spectr a of a slightly higher reso- 
lution than those of Vladilo et al. (1994), resulting in a larger 
number of fitted Na I components. 

In most cases, the Na I column densities from the two opti- 
cal spectra, assigned to our component structure, are compara- 
ble (Table 0). The values disagree for the weak component D^, 
due to different velocity centres of th e fitted components. We 
thus take the value from Vladilo et al. (1994) as an upper limit. 
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A strong discrepancy is visible in the local component Dq. The 
value of Vladilo et al. (1994) is assumed as a lower limit. No 
Nal is detected at +160 kms^^ (component G). 

In the components Db and G, where a column density for 
Zn II cannot be determined, we use instead the linear approxi- 



mation of Vladilo & Centurion ( 1994) 



log 



iV(Znll) 

iV(Nai) 



-0.78(±0.07)logAr(Nai) 



(1) 



-10.31(±0.80), 



based on only Na I. However, these values can only describe 
the general tendency and must not be overinterpreted. 

Temperature T and electron density rig in a gas are related 
to the ionisation ratio of an element X according to 



A^(Xii) ^ Tx 
iV(Xi) a[T) 



(2) 



where Fx is the photoionisation rate and a(T) the recombina- 
tion coefficient. As an initial approximation for Fx, we use the 
values WJl of de Boer et al. (1973 1, based on the mean Galac- 



tic radiation field of Witt & Johnson (1973). Additional contri 



butions fro m cha rge exchange reactions (see e.g. Pequignot & 
Aldrovandi 1986 ) can be neglected, since the density of the gas 
is low. 

Radiative recombination coefficients are taken from 



Pequignot & Aldrovandi ( |198q ). For magnesium we must also 
take into account the dielectronic recombination coefficient 



given by ShuU & van Steenberg (1982), and the contribution 
of dielectronic recombination proces ses at lower temperatures 
calculated by Nussbaumer & Storey ( 1986 ). 

For magnesium we therefore get a temperature dependent 
electron density of 



ne(r) 



FMg iV(Mgi) 
a{T) ' 7V(Mgii)^ 



while we use for sodium 

A/'(Naii) _ Ar(Naii) 7V(Znii) 
iV(Nai) ~ 7V(Znii) ' iV(Nai) 

with 



A^(Naii) _ A^(Na) 
7V(Znll) ~ 7V(Zn) 



50, 



(3) 



(4) 



(5) 



assuming solar abundances. The left approximation in Eq. (g) 
is valid, because in neutral gas the first ionisation stages of 
sodium and zinc are far more abundant than the neutral forms. 

Plotting these two functions into a diagram for each com- 
ponent, we can read the physical conditions in the gas from the 
line curve intersepts (Fig. 0). The resulting limits for tempera- 
tures and electron densities are given in Table 0. Exact limits 
could not be given for all components. 

Only component C yields a well defined intersept that pro- 
vides firm limits on T and rig. For A, a value in the flat part 
of the diagram below 3500 K seems to be most probable. If 
T was significantly higher in this cloud, one would expect a 
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Fig. 2. Electron densities and temperatures of the different ab- 
sorption components. "Mg" denotes the plots calculated from 

7V(Mgi)/7V(Mgii), "Na" those from iV(Nai)/iV(Naii), 
which are derived from A^ (Na i)/N{Zn ii). From the intersepts 
of the curves we estimate the physical parameters of the gas 



much higher state of ionisation and thus much higher electron 
densities in the gas, which would not fit with the diagram. Only 
D;, seems to have a combination of high temperature and high 
electron density. For E and G we probably have temperatures 
between 4000 and 10 000 K. 

Since the curves of magnesium and sodium are nearly par- 
allel below about 3500 K, we need additional information for a 
more precise temperature determination. According to Vladilo 
& Centurion (|1994[), logi?Ar = log[A^(Znll)/A^(Nal)] can be 



used by itself as an indicator for temperature and electron den- 
sity in the gas. In our data the column densities are not exact 
enough for log Rn to provide tight restrictions to these param- 
eters, but we can at least assume cold gas around 100 K for 
log Rn ^0.8 and warm gas for log R^ <; 0.9. 

Under this criterion, we find that component Da and prob- 
ably also F consist of cold gas. While A and C are still un- 
decided, the gas in B is probably warm. For the gas in E and 
especially G we also get a tendency towards higher tempera- 
tures. 

With known electron densities we can now roughly esti- 
mate hydrogen densities and hence the extents of the clouds. 
Assuming the hydrogen in the gas is totally neutral, most of the 
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Table 6. Estimated temperatures, electron densities, hydrogen 
densities, and line of sight extents of the identified velocity 
components 



ne [cm~^] " 



T[K]' 



nu [cm ■*] I [pc] 



A 0.02 - 0.2 

B <0.2 

C 0.08 - 0.32 

Da 0.02 - 32 

Di, < 1.6 

E < 0.06 

F 0.13-4.0 

G < 0.16 



~ 100 - 5000 
3200 - 10000 
5000 - 7000 
~ 100 - 6000 

^5000 
~ 100 - 10000 
~ 100 - 4000 

>3500 



^50'' 2.6 

<?'= > 1.8 

2.7-11'= 10-2 

^50'' ^1.6 

< 50 '^ ;$ 8 

<2'= > 10 

^300* ^0.3 

<5'= ;<0.6 



" from Fig. S 

*" adopting fully neutral H, free electrons mainly from C^ 

" adopting 3% H+ 



electrons stem from carbon, which is the most abundant ionised 
metal in the ISM. For cold gas at about 100 K w e thus c an use 
the carbon abundance of —3.4 dex (de Boer et al. 1987a ) to de- 
rive the hydrogen abundance. This probably is the case in com- 
ponents A, Da, and F. For warm gas of several 1000 K the frac- 
tion of ionised hydrogen can be 1% to 5%, depending on the 
level o f ioni sing radiation in the environment (Vladilo & Cen- 
turion 1994 ; Kulkarni & Heiles 1987 ). Under this assumption, 
electrons from ionised hydrogen are by far dominant. However, 
the wide temperature ranges derived above and the lack of in- 
formation on the radiation field in the clouds make it difficult to 
determine exact densities. For the probably warm components 
B, C, Dfe, E, and G with T ^ 4000 K, we assume a fraction of 
ionised hydrogen of 3%. 

Based on the values derived for nn, we then calculate the 
line of sight extent of each cloud, with appropriate uncertain- 
ties. The results are Usted in Table g[ 

Due to the uncertain ionisation ratio of hydrogen, espe- 
cially in the warm components, the derived densities could be 
smaller than listed in Table ^. McKee & Ostriker (1977) as- 
sume a fractional ionisation of 15% for the warm neutral ISM, 
which would lead to a factor 5 smaller hydrogen densities. In 
the cool components A, Da, and F a non-negligible fraction of 
ionised hydrogen could have an even larger effect on the hy- 
drogen densities. 

As a consequence, also the calculated cloud sizes strongly 
depend on the level of hydrogen ionisation assumed for the 
clouds. The warm clouds are about a factor 5 larger if 15% 
of the hydrogen is ionised. Looking at the cool components, a 
fractional ionisation of only 1 % already means a linear extent 
about a factor 25 larger than stated in Table ra. 



4. Higher ionisation stages towards SN 1993J 



4.1. Column densities 

The gas seen in the absorption lines of higher ions is likely spa- 
tially distinct from the lower ionised gas and has to be treated 
separately. 

For the analysis, we determine apparen t colu mn densities, 
as described by, e.g.. Savage & Sembach ( 1991 ): The optical 
depth T at a wavelength A is approximately 



T(A)~T<,(A)=ln[/o(A)//obs(A)], 



(6) 



with Io{X) the continuum intensity and Jobs (A) the observed 
intensity at this wavelength. For each velocity bin of 20 km s^ ^, 
which is the resolution of the lUE, the logarithmic column den- 
sity \og[Na{v)] apparently occuring in that bin then computes 

as 

log[Naiv)] = logTaiv) - log(/Atr) + 15.876 (7) 

[atomscm~^(20kms"^)"^], 

where / is the transition probability and Ati the wavelength of 
the transition. 

We examine the doublet transitions of Aim (1854 and 
1862 A), Si IV (1393 and 1402 A), and C IV (1548 and 1550 A). 
No absorption of the N V doublet at 1238 and 1242 A is ob- 
served along this line of sight. This lack of N V absorption was 
also observed towards SN 1987A (Fransson et al. 1987) and 



was to be expected. The lower solar abundance of nitrogen and 
the lower /-values for the N V transitions lead to a nominal line 
strength about a factor of 6 lower than that of C IV, making N V 
hard to detect in the halo. Furthermore, the N V transitions are 
located at smaller wavelengths, where the UV flux from the 
supernova had already decreased, thus lowering the signal-to- 
noise ratio in this region of the spectra. 

Rehable absorption profiles are calculated for each of the 
observed spectra. The comparison of the apparent column den- 
sity profiles in each doublet shows no indication of saturation 
for any of the atoms, i.e., no unusual deviation in the peak 
heights. We thus take the mean of the four available profiles 
(two transitions in each of the two spectra) of each ion. Fig. 
g| shows the resulting mean apparent column density profiles. 
C IV showed an unnaturally low column density at —40 km s^^, 
which was disregarded. 

Significantly larger errors for the mean apparent col- 
umn densities at velocities lower than —100 kms^^ 
(Aim: -120 kms^^) and larger than +220 kms^^ (Aim: 
+200 kms^^) indicate a lack of absorption at these veloc- 
ities. Only upper limits can be given here. Especially near 
— 130 kms~^, the rotational velocity of M81 at the position 
of SN 1993 J, no absorption is found. 

We find a broad absorption component around —40 km s^^ 
and another broad absorption between +80 and +200 km s^^, 
the latter probably due to at least 2 clouds of velocities around 
+110 and +170 kms^^. Also, a weak component at kms~^ 
is visible. The integrated column densities for these compo- 
nents are listed in Table ^ 
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Fig. 3. Apparent column densities for the highly ionised atoms 
Civ, Si IV, and Aim. The extremely low Civ value at 
—40 km s^^ was disregarded. Where no column density is plot- 
ted, the profile indicated upper limits lower than the end points 
of the profiles 

Table 7. High ion metal column densities towards SN 1993 J. 
Only upper limits can be given for gas at —130 kms^^. The 
errors for the other values are in the range of 0.15 dex. In the 
rightmost column the total column density for absorption above 
+80kms"i is listed 



log [J2^ Niv)] [cm-2] 


i;LSR[kms-^] -130 


-40 +110 +170 


110 & 170 


CIV < 13.1 
Si IV < 12.8 
Aim < 12.3 


14.3 ~ 13.0 14.0 13.9 

14.0 ~ 12.8 13.7 13.7 

13.1 - 12.9 12.9 


14.3 
14.0 
13.2 



4.2. Physical conditions in the gas 

lonisation of gas to the level of Si IV and Al III can be explained 
by stellar photons. For C IV the necessary ionisation energy is 
too large to be provided by stellar radiation, but these ions can 
be easily produced, for example, in the cooling process of a 
galactic fountain. 

Looking at the apparent column density profiles, we find 
similar profiles for Si IV and C IV in the negative velocity range, 
with a nearly constant ratio log[iV(Civ)/N(Silv)] of about 
2.5. The profile of Al III differs, but this might be due to the 
weakness of the absorption. Both C IV and Si IV possibly ex- 
ist in the same spatial environment, in a border region between 
hot, infalling fountain gas and cooler, photoionised or recom- 
bined, gas. This spatial coexistence seems to be a common phe- 



nomenon in halo gas, it was also found by Savage et al. (1989) 
on the line of sight to SN 1987 A. 

Nothing comparable is found in the positive velocity range. 
The profile of C IV clearly differs from that of Si IV. We proba- 
bly have spatially distinct clouds at similar velocities here, pos- 
sibly in the vicinity of M 8 1 . 



5. Spatial arrangement of the clouds 

An estimate of the locations of the absorbing clouds is difficult 
for the highly ionised gas. However, we can distinguish three 
basic groups of gas; gas within the Milky Way, either in the 
disk or in its halo, gas in M81, and intergalactic gas at some 
position in the space between M 8 1 and the Galaxy, probably 
near the M 81 group. 

5.7. Gas in the Milky Way 

In the low ion transitions, components C, Dq, and D^ are con- 
sidered to be Galactic gas. Their velocities and their clear vis- 
ibility in the 21 cm radio emission spectra let any other inter- 
pretation than nearby gas seem improbable. 

Component C, from its velocity of —50 km s^^, is by defi- 
nition an IVC in our Galaxy, probably part of the L ow La titude 
Intermediate Velocity (LLIV) Arch (Kuntz & Danly |l996| ). The 
large 6-value and column density indicate a possible blend of 
more than one cloud. 

Component D^ shows the typical parameters for cold local 
disk gas. 

Component Df, may represent absorption by gas of the Lo- 
cal Cloud, in which our solar system is imbedded. F or Lo - 
cal Cloud gas in the direction of M81, Genova et al. ( 199C| ) 
find a temperature around 11500 K and a hydrogen density of 
uh — 0.1 cm^'^ at a velocity of +11 kms~^. The magnesium 
column density has a low value of only log A'^(Mg) < 12.1. 
The extent of this gas is about 3 — 4 pc. These parameters agree 
nicely with our results. 

In the profiles of the highly ionised atoms we find ad- 
ditional Milky Way absorption components: The absorption 
around kms^^ has its origin probably in Galactic disk gas, 
possibly in gas of the Local Bubble, surrounding the Local 
Cloud. The gas seen around —40 kms~^ could be related to 
the IVC component C in neutral gas. A transition between neu- 
tral and ionised gas also would explain the coexistence of C IV 
and Si IV found in Chapter i.2 . We obviously have no contribu- 
tion from the nearby Galactic HVC Complex C, which, in this 
general direction, shows velocities of — 100 kms^^ and below 
(Wakker|l99l||. 



5.2. Gas in M 81 

We find neutral gas in absorption at velocities near 
— 130 kms^^. From the radiointerferometry map of M81 gas 



by Rots (1975), this is the rotational velocity of the galaxy at 



the position of the supernova. Peak-shaped emission is not visi- 
ble in the H I 21 cm emission spectra at this velocity, but rather 
only broad background emission from the disk of M 8 1 . The 
gas in absorption is thus probably located at a large distance, at 
which it could not be detected separately in the 9' radio beam 
of the Effelsberg telescope. For these reasons the components 
A and B can be assumed to be located in the disk or the lower 
haloofM81. 
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In component A, the H l column density is comparable to 
local gas. The metal column densities in this component are 
very similar to those in Milky Way gas. We probably see cold 
M81 disk gas at -130 km s"^ here. 

Also assuming a location in M81, component B can be 
interpreted as an IVC or a blend of such clouds, at a velocity of 
- 40 km s " ^ relative to the M 8 1 disk. 

No highly ionised gas is detected that would have to be 
associated with a location near M 81. 



M81 




5. 



,5. Intergalactic gas in front ofMSl 



The three absorption components E, F, and G seen at positive 
velocities can hardly be assigned to the Galaxy or M 81 itself. 
They must be intergalactic clouds, located probably in the IGM 
of the M 8 1 group. The gas seen in the highly ionised atoms 
at positive velocities is associated with intergalactic clouds as 
well. 

The hydrogen column density of the cooler component E 
is comparable to that of the IVC components B and C, with 
metal depletions in a normal range. 

In the cold component F we find a hydrogen column den- 
sity about one order of magnitude smaller than for the local 
Galactic component Dq. The depletion is low, with presumably 
no depletion at all for oxygen. This is an indicator for a low 
dust content in the gas. If we assume a location in the vicinity 
of M81, the interpretation as an infalling M81 HVC is possi- 
ble. The low temperature, however, contradicts an origin from 
a galactic fountain mechanism. 

Component G shows warm gas of very low hydrogen col- 
umn density and low metal depletion. Similar to component E, 
it could be high velocity gas still in its cooling phase, infalling 
towards M81. 

We find very low electron densities for the warm neutral 
gas components E and G at positive velocities. This is unusual, 
since our calculations are based on the Galactic disk radiation 
field, which is too large for halo gas. The disk radiation fields 
of the Galaxy and M 8 1 should be comparable. However, if the 
components are located at some distance from the disk of M 8 1 , 
as assumed from their velocities, the level of radiation should 
be significantly lower, leading to even lower electron densities 
than derived above. We thus conclude that an enhanced level 
of ionising radiation exists at the location of the components 
E and G, at least at the level of a disk radiation field. Evi- 
dence for enhanced ionising radiation in a range typical for the 
Galactic disk was also inferred by Vladilo et al. (1994) from 
the iV(Caii)/iV(Cai) ratio. 

Radiation from hot disk stars is an improbable source of 
energy for this level of ionisation. The interpretation of typical 
HVCs seems unlikely under these conditions. If located in the 
foreground of M 81, the dynamic intergalactic medium of the 
M 8 1 group itself could provide the energy needed. The source 
of energy is unclear, but could be shocks or plasma effects like 
magnetic reconnection (Birk et al. 1998). Also mixing effects 



MW 
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Fig. 4. Schematical view of the lines of sight from the Milky 
Way to the two nearby supernovae. In the lower right the Mag- 
ellanic Clouds are visible at about 50 kpc distance; M81 at 
3.63 Mpc is found in the upper left corner 

at the border region between hot intergalactic gas and cooler 
tidal debris are possible. 

The low pressure in the IGM would explain the probably 
large extent of the hot component E. Component F, from its 
velocity assigned to this intergalactic gas, is possibly the cold, 
shielded core of some intergalactic cloud. This would explain 
the high hydrogen density at a low extent. The high density 
nn > 300 cm~'^ derived assuming free electrons solely from 
carbon might indicate that at least a small fraction of hydrogen 
is ionised in this component. This would mean a lower density 
and a larger extent, rather typical for an extended cold cloud. 

The gas seen in the high ion absorption lines at velocities 
larger than +80 km s^^ must also be located in front of M 8 1 . 



Bregman ( 1980 ) predicts mainly negative velocities in Galac- 
tic neutral gas in the direction of M 8 1 . Of course, the highly 
ionised gas cannot be identical to this neutral gas. Nevertheless, 
if both ionisation stages are only different states of infalling 
gas, their velocities should be comparable. Positive velocities 
around +140 kms^^, however, have been found by Yun et al. 



( |1993D in the tidal tail region connecting M 81 and M 82. 

Since we find indicators for a high energy environment 
from the lower ionisation stages, the coexistence of highly and 
low ionised gas and cold, shielded cloud cores in the M81 
group IGM is plau sible. S imilar conditions were recently found 
by Sembach et al. (2000) in the Magellanic Stream and in outer 
Galaxy HVCs, where neutral H I gas is visible together with 
O VI at the interface to the surrounding hot medium. 

6. Combined SN 1987A to SN 1993J line of sight 

6.1. Line of sight and velocity structure 

The antipodal directions of the two supernovae allow us to 
combine their lines of sight into a single one that traverses three 
galaxies in a direction nearly perpendicular to the plane o f the 
Local Group (Fig. |). As noted by de Boer et al. ( 1987b| ), the 
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Fig. 5. Absorption profiles created for the combined line of 
sight from SN 1987 A to SN 1993 J, based on the spectra to each 
of the supernovae. SN 1987A to SN 1993J. The velocity origin 
is set to the position of SN 1987 A. The marks in the top of 
the figure indicate which absorption structure is due to which 
of the individual gas components along the total line of sight. 
They now heavily blend, of course 



position of SN 1987A is on the backside of the LMC. Starting 
there, the line of sight crosses the local LMC gas, the HVCs and 
IVCs of the southern Galactic halo, the full disk of the Milky 
Way, and its northern halo with the IVC. It continues on, sam- 
pling the intergalactic medium towards the M 8 1 group and the 
halo of M 8 1 , before it reaches the disk of M 8 1 itself at the lo- 
cation of SN 1 993 J. The total length of this line of sight is about 
3.7 Mpc, sampling gas at a total velocity width of 500 km s^^. 
To create a spectrum over the combined line of sight we 
merge the absorption line spectra of the two supernovae. We 
examine the transitions of Sill at 1304, 1526 and 1808 A, as 
well as the C IV transition at 1548 A. For the combination we 
use the ripple corrected spectra and sum over the absorption in 
each velocity bin. We do lose some equivalent width in those 
absorption components which are saturated. It is shown below 
(Sect. PjA}, however, that this loss has a negligible effect on 
the total equivalent width. It is not neccessary for this exami- 
nation to add the optical depths and generate Voigt profiles of 
the absorption, which is the mathematically correct procedure. 
Fig. shows the combined absorption profiles for the total line 
of sight. We see a complex blending of absorptions from the 
different galaxies in the absorption profile. Without knowledge 
about the individual galaxies it is impossible to associate indi- 
vidual absorption components with a particular galaxy. 



6.2. Column densities 

We obtain the equivalent widths for the complete absorption 
structure by a trapezium fit. For Si II we can fit these values 
to a theoretical curve of growth with a doppler parameter of 
6 ~ 80 kms^^. A 6-value that large was to be expected for 
such massive blending of absorption components. The total ab- 
sorption leads to a column density of log Aftotai(Si ii) = 16.1. 
If we just sum over the equivalent widths of the single compo- 
nents and fit these to a curve of growth, we get the same value 
of logiVtotai(Siii) == 16.1, but with a 6-value of 100 kms~^ 
Obviously, the column density for these data is not very sensi- 
tive to changes in the equivalent width. This justifies the proce- 
dure of just adding the absorption profiles instead of correctly 
adding the optical depths. 

From the sum over the single column densities we get 
log A'totai(Siii) = 16.6, which is about a factor of 3 larger 
than from the curve of growth. This deviation shows that, in 
absorption profiles with heavy line blending, the total column 
density can easily be underestimated by a significant amount. 

No reliable total column density can be determined for 
C IV. With only one data point we cannot use the curve of 
growth method. The apparent column density method is not 
suitable here, due to the saturation of the profile. The sum over 
the single column densities is lo g A^tot ai(Civ) ^ 15.0, using 
the value of 14.7 (Savage et al. |1989| ) for SN 1987A and the 
values from Table || for SN 1993J. 

7. Comparison with quasar lines of sight 

Complex absorption structures seen in the spectra of quasars 
are commonly identified with absorption by galaxies, galactic 
halos, or other diffuse objects such as dwarf galaxies along the 



line of sight (s ee, e.g .. Savage & Jeske 1981, York et al. 1986, 
Yanny & York 1992 ). We now compare typical quasar absorp- 
tion with the absorption structure found on our Local Group 
line of sight. We use equivalent widths listed for three different 
quasar sight lines from Savage & Jeske (1981 1 and those of on e 
quasar in the Hubble Deep Field South from Savaglio ( 1998 ). 
The velocity widths of the absorption complexes are between 
150 and 550 kms^^, which is comparable to the 500 kms^^ 
found for the local absorption. The data are listed in Table |[ 

The column densities in Sill are generally lower on the 
quasar sight lines than in the SN 1987 A - SN 1993 J gas, while 
in C IV the values are similar The relation of lower to higher 
ionisation stages is thus significantly lower on the quasar sight 
lines. Both results are expected, since the local line of sight is 
an extreme case of three galactic disks along one line of sight. 
Most of the neutral gas is found in the disks of galaxies, so this 
explains the high column density in Si II. The halos of galax- 
ies are far more extended than the disks, so the probability of 
sampling halo gas in a quasar sight line is far higher than for 
disk gas. This leads to the observed relation of lower to higher 
ionisation stages in the quasar absorption complexes and to the 
generally lower column density in the less dense halo gas. 

Dwarf galaxy gas or gas of any other diffuse nature may 
also produce absorption lines. Substantial Si II may be present 
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Table 8. Quasar column densities in comparison to column 
densities on the line of sight from SN 1987 A to SN 1993 J 



log iV [cm~^] 



Note' 







Sill 


Civ 


Civ 
Sin 




QSO 1756+237 


1.673 


> 14.3 


^ 15.8 


< 1.5 


a 


QSO 0002-422 


2.302 


14.9 


14.9 


0.0 


b 


PKS 2126-158 


2.638 


14.8 


> 15.0 


>0.2 


c 




2.769 


14.3 


14.7 


0.4 


c 


QSOJ2233-606 


1.787 


< 14.9 


14.3 


> -0.6 


d 




1.869 


< 14.8 


14.5 


> -0.3 


d 




1.928 


< 13.6 


< 13.7 


~0.1 


d 




1.943 


< 14.0 


< 14.8 


^ 0.8 


d 




2.077 


< 13.6 


13.3 


> -0.3 


d 




2.198 


< 13.3 


13.8 


> 0.5 


d 




2.206 


< 14.3 


14.2 


> -0.1 


d 



SN 1987A to SN 1993J 


> 16.1 


;> 15.0 


<. 


-1.1 


LMC gas 


16.4* 


> 14.8'= 


> 


-1.6 


Clouds F+G to SN 1993 J 


14.5 


14.3 




-0.2 



Data from (a) Tu rnshek et al. ( 197S )\(h) Sargent et al. ( 1979 1; (c) 



Young et al. ( |1979| ); (d) Savaglio ( |1998| ) 

'' Gas at 7JLSR > 190 km s~\ from Welty et al. ( |l999|) 

" Gas at wlsr > 190 km s-\ from Savage et al.^l989|) 



in dwarf galaxies, as in the case of the LMC. Free gas outside 
galaxies containing metals is most likely material expelled by 
such galaxies. 

The ratio A^(Clv)/A^(Sill) could possibly be used as an 
indicator for halo gas along lines of sight. For local gas we get 
a value of < —1.1 dex, while the values for quasar absorp- 
tion complexes range between —0.6 and 1.5 dex. This ratio can 
be influenced by different depletions of the elements, however, 
especially since both ionisation stages exist in completely dif- 
ferent environments. 

In particular, we can compare the cloud complex Fh-G lo- 
cated in the intergalactic medium in front of M 8 1 with these 
results. Since we do not detect high ions at the velocity of M 8 1 
component E, we consider only the ratio for components F and 
G alone (see Table 0). In contrast to the total line of sight, the 
ratio of N{C lv)/A^(Si ll) for the intergalactic cloud complex 
is very similar to those from QSO absorption systems, with 
comparably low column densities. 

8. The intergalactic clouds in front of M 81 

The cloud complex Fh-G found in the foreground of M 8 1 does 
not show a low density, metal poor structure as one usually 
assumes for intergalactic gas. Instead, it shows a mixture of 
cold and warm, normal density, metal rich gas, in coexistence 
with highly ionised hot components, as normally expected on 
lines of sight crossing a galactic disk. This indicates that the 
Fh-G complex probably contains gas of galactic disk origin. 

The complex extragalactic gas features of the M 8 1 group 
have been mapped in emission of both Hi (see, e.g., Cottrell 



1977| ; Appleton et al. |1981| ; Apple ton & van der Hulst |1988| ) 
and CO (e.g., Brouillet et al. 1991 ). Identifying gas as lying in 



front of the bright M 8 1 disk, however, is difficult from emis- 
sion. The existence of a separate cloud would follow from the 
difference of its velocity with respect to the velocity of the gas 
of the background object (here M 81), either in H I or in CO. 
The UV absorption lines show that the detected gas is indeed 
in front of M 8 1 . Furthermore, they provide us with full infor- 
mation on metallicities and ionisation ratios in the intergalactic 
gas in this foreground region. 

Molecular gas was found in the intergalactic medium sur- 
rounding M81 (Brouillet et al. 



1992; Walter & Heithausen 



1999). We do not know if molecules exist in component F+G, 



since no radio molecules were looked for and the visual spec- 
tra of SN 1993J gave only upper limits for CH+ (Vladilo et al. 



1994). Given the low temperature of absorption component F 



(see Sect. 3.2), the existence of molecular gas in the F+G cloud 



complex is possible. 

The origin of the gas on our line of sight is possibly M 82. 
The absorption velocities are comparable to the ones found by 
Yun et al. ( |1993[ ) for the Hi tidal tail co nnecting M81 and 
M 82. A recent simulation by Sofue (119981) indicates that M 82 



could have lost its complete disk in an encounter with M 81. If 
this were so, we would expect processed disk gas in the fore- 
ground of M 81, possibly also containing stars. The Sofue sce- 
nario could thus explain the metallicities as well as the galaxy 
disk-like ionisation structure of the gas. We note here, that an- 
other stripping event may ha ve taken place in the M 8 1 group. 
Walter & Heithausen (1999) found CO gas displaced with re- 
spect to the optically well-defined galaxy NGC 3077, i.e., here 
without indications for associated stars. 

Our procedure to derive electron densities in the gas com- 
ponents is based on the average Milky Way disk radiation field. 
The level of radiation should be different in intergalactic space 
far away from galaxy disks. Nevertheless, the intergalactic ab- 
sorption components show signs for ionisation comparable to 
that in the radiation field of the Galaxy or even larger. Another 
explanation for the observed level of ionisation could be the 
presence of a stellar component in the gas, not unrealistic in a 
scenario of material ripped from M 82 in a close encounter. If 
M 82 indeed lost its spiral arms, stars could well be present in 
the intergalactic medium of the M 8 1 group. 

The detection of highly ionised gas in the complex F+G 
demonstrates that neutral and ionised gas are mixed, very simi- 
lar to what is seen in QSO absorption line systems. It is unclear, 
however, just where exactly the ionised gas is situated. It is ei- 
ther in the border region of the hot intergalactic plasma, ther- 
mally excited by collisions, or it is gas excited by a radiation 
source somewhere within the M 8 1 system. 
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